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We discuss exact analytical solutions of a variety of statistical models recently obtained for finite systems by a novel powerful mathemati- 
cal method, the Laplace-Fourier transform. Among them are a constrained version of the statistical multifragmentation model, the Gas 
of Bags Model and the Hills and Dales Model of surface partition. Thus, the Laplace-Fourier transform allows one to study the nuclear 
matter equation of state, the equation of state of hadronic and quark gluon matter and surface partitions on the same footing. A 
complete analysis of the isobaric partition singularities of these models is done for finite systems. The developed formalism allows us, for 
the first time, to exactly define the finite volume analogs of gaseous, liquid and mixed phases of these models from the first principles of 
statistical mechanics and demonstrate the pitfalls of earlier works. The found solutions may be used for building up a new theoretical 
apparatus to rigorously study phase transitions in finite systems. The strategic directions of future research opened by these exact 
results are also discussed. 



1. Theoretical Description of Phase Transitions in Finite Systems 

A rigorous theory of critical phenomena in finite systems was not built up to now. However, the experimental 
studies of phase transitions (PTs) in some systems demand the formulation of such a theory. In particular, the 
investigations of the nuclear liquid-gas PT [iHS] require the development of theoretical approaches which would 
allow us to study the critical phenomena without going into the thermodynamic limit — > ex: (1^ is the volume of 
the system) because such a limit does not exist due the long range Coulomb interaction. Therefore, there is a great 
need in the theoretical approaches which may shed light on the "internal mechanism" of how the PTs happen in 
finite systems. 

The general situation in the theory of critical phenomena for finite (small) systems is not very optimistic at the 
moment because theoretical progress in this field has been slow. It is well known that the mathematical theory of 
phase transitions was worked out by T. D. Lee and C. N. Yang j^. Unfortunately, there is no direct generic relation 
between the physical observables and zeros of the grand canonical partition in a complex fugacity plane. Therefore, 
we know very well what are the gaseous phase and liquid at infinite volumes: mixture of fragments of all sizes and 
ocean, respectively. This is known both for pure phases and for their mixture, but, despite some limited success [S], 
this general approach is not useful for the specific problems of critical phenomena in finite systems (see Sect. VIII 
below) . 

The tremendous complexity of critical phenomena in finite systems prevented their systematic and rigorous 
theoretical study. For instance, even the best formulation of the statistical mechanics and thermodynamics of finite 
systems by Hill [6] is not rigorous while discussing PTs. As a result, the absence of a well established definition of 
the liquid and mixed phase for finite volumes delays the progress of several related fields, including the theoretical 
and experimental searches for the reliable signals of several PTs which are expected to exist in strongly interacting 
matter. Therefore, the task of highest priority of the theory of critical phenomena is to define the finite volume 
analogs of phases from first principles of statistical mechanics. At present it is unclear whether such definitions can 
be made for a general case, but it turns out that such finite volume definitions can be formulated for a variety of 
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realistic nonclassical (= non mean-field) statistical models which are successfully used in nuclear multifragmentation 
and in relativistic heavy coUsisions. 

About 25 years ago, when the theoretical foundations of nuclear multifragmentation were established, there was 
an illusion that the theoretical basis is simple and clear and, therefore, we need only the data and models which 
will describe them. The analysis of finite volume systems has proven to be very difficult. However, there was a clear 
way out of troubles by making numerical codes that are able to describe the data. This is, of course, a common way 
to handle such problems and there were many successes achieved in this way [T}]31[7j. However, there is another side 
of the coin which tells us that our understanding did not change much since then. This is so because the numerical 
simulations of this level do not provide us with any proof. At best they just demonstrate something. With time 
the number of codes increased, but the common theoretical approach was not developed. This led to a bitter result 
- there are many good guesses in the nuclear multifragmentation community, but, unfortunately, little analytical 
work to back up these expectations. As a result the absence of a firm theoretical ground led to formulation of such 
highly speculative "signals" of the nuclear liquid- vapor PT as negative heat capacity [8, 9J, bimodality [lOj, which 
later on were disproved, in Refs flTj and [12], respectively. 

Thus, there is a paradoxic situation: there are many experimental data and facts, but there is no a single 
theoretical approach which is able to describe them. Similar to the searches for quark-gluon plasma (QGP) |13| 
there is lack of a firm and rigorous theoretical approach to describe phase transitions in finite systems. 

However, our understanding of the multifragmentation phenomenon [hri^ was improved recently, when an exact 
analytical solution of a simplified version of the statistical multifragmentation model (SMM) jl4l[T5] was found in 
Refs. [lOI'lT]. These analytical results not only allowed us to understand the important role of the Fisher exponent 
T on the phase structure of the nuclear liquid-gas PT and the properties of its (tri)critical point, but to calculate 
the critical indices a',/3,7',5 of the SMM [18] as functions of index r. The determination of the simplified SMM 
exponents allowed us to show explicitly [18] that, in contrast to expectations, the scaling relations for critical indices 
of the SMM differ from the corresponding relations of a well known Fisher droplet model (FDM) 119] . This exact 
analytical solution allowed us to predict a narrow range of values, 1.799 < t < 1.846, which, in contrast to FDM 
value TpDM ~ 2.16, is consistent with ISiS Collaboration data [50] and EOS Collaboration data [5T]. This finding 
is not only of a principal theoretical importance, since it allows one to find out the universality class of the nuclear 
liquid-gas phase transition, if r index can be determined from experimental mass distribution of fragments, but also 
it enhanced a great activity in extracting the value of r exponent from the data [22J . 

It is necessary to stress that such results in principle cannot be obtained either within the widely used mean- 
filed approach or numerically. This is the reason why exactly solvable models with phase transitions play a special 
role in statistical mechanics - they are the benchmarks of our understanding of critical phenomena that occur in 
more complicated substances. They are our theoretical laboratories, where we can study the most fundamental 
problems of critical phenomena which cannot be studied elsewhere. Their great advantage compared to other 
methods is that they provide us with the information obtained directly from the first principles of statistical 
mechanics being unspoiled by mean-field or other simplifying approximations without which the analytical analysis 
is usually impossible. On the other hand an exact analytical solution gives the physical picture of PT, which cannot 
be obtained by numerical evaluation. Therefore, one can expect that an extension of the exact analytical solutions 
to finite systems may provide us with the ultimate and reliable experimental signals of the nuclear liquid-vapor PT 
which are established on a firm theoretical ground of statistical mechanics. This, however, is a very difficult general 
task of the critical phenomena theory in finite systems. 

Fortunately, we do not need to solve this very general task, but to find its solution for a specific problem of 
nuclear liquid-gas PT, which is less complicated and more realistic. In this case the straightforward way is to 
start from a few statistical models, like FDM and/or SMM, which are successful in describing the most of the 
experimental data. A systematic study of the various modifications of the FDM for finite volumes was performed by 
Moretto and collaborators [23] and it led to a discovery of thermal reducibility of the fragment charge spectra |3], 
to a determination of a quantitative liquid-vapor phase diagram containing the coexistence line up to critical 
temperature for small systems [24, 25j , to the generalization of the FDM for finite systems and to a formulation 
of the complement concept [5S1[27] which allows one to account for finite size effects of (small) liquid drop on the 
properties of its vapor. However, such a systematic analysis for the SMM was not possible until recently, when its 
finite volume analytical solution was found in ^28j. 
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An invention of a new powerful mathematical method |28) . the Laplace- Fourier transform, is a major theoretical 
breakthrough in the statistical mechanics of finite systems of the last decade because it allowed us to solve exactly 
not only the simplified SMM for finite volumes [28J, but also a variety of statistical surface partitions for finite 
clusters [53] and to find out their surface entropy and to shed light on a source of the Fisher exponent r. It was 
shown [28] that for finite volumes the analysis of the grand canonical partition (GCP) of the simplified SMM is 
reduced to the analysis of the simple poles of the corresponding isobaric partition, obtained as a Laplace-Fourier 
transform of the GCP. Such a representation of the GCP allows one not only to show from first principles that 
for finite systems there exist the complex values of the effective chemical potential, but to define the finite volume 
analogs of phases straightforwardly. Moreover, this method allows one to include into consideration all complicated 
features of the interaction (including the Coulomb one) which have been neglected in the simplified SMM because 
it was originally formulated for infinite nuclear matter. Consequently, the Laplace-Fourier transform method opens 
a principally new possibility to study the nuclear liquid-gas phase transition directly from the partition of finite 
system without taking its thermodynamic limit. Now this method is also applied [30] to the finite volume formulation 
of the Gas of Bags Model (GBM) [31] which is used to describe the PT between the hadronic matter and QGP. 
Thus, the Laplace-Fourier transform method not only gives an analytical solution for a variety of statistical models 
with PTs in finite volumes, but provides us with a common framework for several critical phenomena in strongly 
interacting matter. Therefore, it turns out that further applications and developments of this method are very 
promising and important not only for the nuclear multifragmentation community, but for several communities 
studying PTs in finite systems because this method may provide them with the firm theoretical foundations and a 
common theoretical language. 

It is necessary to remember that further progress of this approach and its extension to other communities 
cannot be successfully achieved without new theoretical ideas about formalism it-self and its applications to the 
data measured in low and high energy nuclear collisions. Both of these require essential and coherent efforts of two 
or three theoretical groups working on the theory of PTs in finite systems, which, according to our best knowledge, 
do not exist at the moment either in multifragmentation community or elsewhere. Therefore, the second task of 
highest priority is to attract young and promising theoretical students to these theoretical problems and create 
the necessary manpower to solve the up coming problems. Otherwise the negative consequences of a complete 
dominance of experimental groups and numerical codes will never be overcome and a good chance to build up a 
common theoretical apparatus for a few PTs will be lost forever. If this will be the case, then an essential part of 
the nuclear physics associated with nuclear multifragmentation will have no chance to survive in the next years. 

Therefore, the first necessary step to resolve these two tasks of highest priority is to formulate the up to 
day achievements of the exactly solvable models and to discuss the strategy for their further developments and 
improvements along with their possible impact on transport and hydrodynamic approaches. For these reasons the 
paper is organized as follows: in Sect. II we formulate the simplified SMM and present its analytical solution in 
thermodynamic limit; in Sect. Ill we discuss the necessary conditions for PT of given order and their relation to 
the singularities of the isobaric partition and apply these findings to the simplified SMM; Sect. IV is devoted to 
the SMM critical indices as the functions of Fisher exponent r and their scaling relations; the Laplace-Fourier 
transform method is presented in Sect. V along with an exact analytical solution of the simplified SMM which 
has a constraint on the size of largest fragment, whereas the analysis of its isobaric partition singularities and the 
meaning of the complex values of free energy are given in Sect. VI; Sect. VII and VIII are devoted to the discussion 
of the case without PT and with it, respectively; at the end of Sect. VIII there is a discussion of the Chomaz and 
Gulminelli's approach to bimodality [5J; in Sect. IX we discuss the finite volume modifications of the Gas of Bags, 

1. e. the statistical model describing the PT between hadrons and QGP, whereas in Sect. X we formulate the Hills 
and Dales Model for the surface partition and present the limit of the vanishing amplitudes of deformations; and, 
finally, in Sect. XI we discuss the strategy of future research which is necessary to build up a truly microscopic 
kinetics of phase transitions in finite systems. 

2. Statistical Multifragmentation in Thermodynamic Limit 

The system states in the SMM are specified by the multiplicity sets {n^} (nj. = 0, 1, 2, ...) of fc-nucleon fragments. 
The partition function of a single fragment with k nucleons is [ll: V(f>kiT) = V (mTfc/27r)'^^^ Zk , where k = 1,2, A 
(A is the total number of nucleons in the system), V and T are, respectively, the volume and the temperature of the 
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system, m is the nucleoli mass. The first two factors on the right hand side (r.h.s.) of the single fragment partition 
originate from the non-relativistic thermal motion and the last factor, Zk, represents the intrinsic partition function 
of the fc-nucleon fragment. Therefore, the function 4>k{T) is a phase space density of the k-nucleon fragment. For 
k — I (nucleon) we take zi = 4 (4 internal spin-isospin states) and for fragments with fc > 1 we use the expression 
motivated by the liquid drop model (see details in Ref. [IJ): Zk = exp(— /fe/T), with fragment free energy 

fk = -WiT)k + a{T) + (r + 3/2)T In k , (1) 

with W{T) — Wo + T^/eq. Here Wo — 16 MeV is the bulk binding energy per nucleon. T^/eo is the contribution 
of the excited states taken in the Fermi-gas approximation (eq — 16 MeV). (j(r) is the temperature dependent 
surface tension parameterized in the following relation: cr(T) — a{T)\sMM = o'o[(T)? — T^)/{T^ + T^)]^/^, with 
o-o = 18 MeV and Tc = 18 MeV (a ^ at T > Tc). The last contribution in Eq. Q involves the famous Fisher's 
term with dimensionless parameter r. As we will show later, at the critical (tricritical) point the fragment mass 
distribution will lose it exponential form and will become a power law fc""^. 

It is necessary to stress that the SMM parametrization of the surface tension coefficient is not a unique one. For 
instance, the FDM successfully employs another one cr(T)|FDM — '^o[^ ~ T/Tc]. As we shall see in Sect. IV the 
temperature dependence of the surface tension coefficient in the vicinity of the critical point will define the critical 
indices of the model, but the following mathematical analysis of the SMM is general and is valid for an arbitrary 
cr(T) function. 

The canonical partition function (CPF) of nuclear fragments in the SMM has the following form: 



{"fc} 



k=l 



S{A - kuk) . (2) 



In Eq. ^ the nuclear fragments are treated as point-like objects. However, these fragments have non-zero proper 
volumes and they should not overlap in the coordinate space. In the excluded volume (Van der Waals) approximation 
this is achieved by substituting the total volume V in Eq. ([2]) by the free (available) volume Vf = V — kuk, 
where b = 1/po {po = 0.16 fm~'^ is the normal nuclear density). Therefore, the corrected CPF becomes: Za{V,T) = 
Zjf (F - bA,T). The SMM defined by Eq. ^ was studied numerically in Refs. [T4l[T5] . This is a simplified version 
of the SMM, since the symmetry and Coulomb contributions are neglected. However, its investigation appears to 
be of principal importance for studies of the nuclear liquid-gas phase transition. 

The calculation of Za{V,T) is difficult due to the constraint J2k ^"-fc ~ ^- This difficulty can be partly avoided 
by evaluating the grand canonical partition (GCP) 

oo 

Z{V, T, m) = ^ exp (^) Za{V, T) Q{V - hA) , (3) 

where p denotes a chemical potential. The calculation of Z is still rather difficult. The summation over {n^} sets in 
Za cannot be performed analytically because of additional ^-dependence in the free volume Vf and the restriction 
Vf > 0. The presence of the theta- function in the GCP ([3]) guarantees that only configurations with positive value 
of the free volume are counted. However, similarly to the delta function restriction in Eq. ([2]), it makes again the 
calculation of Z{V,T,^) ^ to be rather difficult. This problem was resolved [inilTT] by performing the Laplace 
transformation of Z{V,T,^). This introduces the so-called isobaric partition function (IP) [31]: 

poo pQC 1 \ 7lL. 

"'0 "'0 {„fc} k 

dV e--^' exp \v'Y,<\>k e^^^ \ . (4) 



k=\ 



After changing the integration variable V ^ V' , the constraint of 0-function has disappeared. Then all Uk were 
summed independently leading to the exponential function. Now the integration over V' in Eq. ^ can be done 
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resulting in 
where 



T{s,T,n) = exp 



(fi — sbT)k 



k=l 



T 



fmTY 



zi exp 



sbT 



T 



exp 



fe=2 



{fi - shT)k - ak^'^ 



T 



(5) 

(6) 



Here we have introduced the shifted chemical potential jl = ^ + W{T). From the definition of pressure in the 
grand canonical ensemble it follows that, in the thermodynamic limit, the GCP of the system behaves as 



T lim 



In Z(F,r,/x) 
V 







^ exp 

y^oo 


T 



(7) 



An exponentially over V increasing part of Z{V^T,fi) in the right-hand side of Eq. ([7]) generates the rightmost 
singularity s* of the function Z{s,T,fi), because for s < p{T,fi)/T the ^-integral for Z{s,T,^) ^ diverges at 
its upper limit. Therefore, in the thermodynamic limit, V ^ oo the system pressure is defined by this rightmost 
singularity, s*(T,^), of IP Z{s,T,fj.) (g]): 



piT,fi) = Ts*{T,ti) 



(8) 



Note that this simple connection of the rightmost s-singularity of Z, Eq. (j4]), to the asymptotic, V ^ oo, behavior 
of Z, Eq. (O; is a general mathematical property of the Laplace transform. Due to this property the study of the 
system behavior in the thermodynamic limit V ^ oo can be reduced to the investigation of the singularities of Z. 

3. Singularities of Isobaric Partition and Phase Transitions 

The IP, Eq. has two types of singularities: 1) the simple pole singularity defined by the equation 

Sg{T,n) ^ T{sg,T,n) , (9) 

2) the singularity of the function T{s, T, fj.) it-self at the point si where the coefficient in linear over k terms in the 
exponent is equal to zero. 



si{T,fi) = ^ . (10) 
The simple pole singularity corresponds to the gaseous phase where pressure is determined by the equation 



/ mT 



3/2 



T 



zi exp 



M - bpg 
T 



exp 



k=2 



[fi - bpg)k - CTfc2/3 



T 



The singularity si{T, /i) of the function J^(s, T, /i) 1^ defines the liquid pressure 

Pi{T,n) = T si{T,fi) 



b 



(11) 



(12) 



In the considered model the liquid phase is represented by an infinite fragment, i.e. it corresponds to the 
macroscopic population of the single mode k = oo. Here one can see the analogy with the Bose condensation where 
the macroscopic population of a single mode occurs in the momentum space. 

In the (T, yu)-regions where fl < bpg{T, /i) the gas phase dominates (pg > pi), while the liquid phase corresponds 
to /i > bpg{T,fj,). The liquid-gas phase transition occurs when two singularities coincide, i.e. Sg{T,fi) = si{T,fi). A 
schematic view of singular points is shown in Fig. la for T < Tc, i.e. when a > 0. The two-phase coexistence region 
is therefore defined by the equation 



PiiT,fi) = Pg{T,n) , i.e., jl = bpg{T,fi) 



(13) 
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One can easily see that J-{s,T,fi) is monotonously decreasing function of s. The necessary condition for the phase 
transition is that this function remains finite in its singular point s; — fl/Tb: 



T{si,T,fj.) < oo 



(14) 



The convergence of T is determined by t and tr. At r = the condition (fT4|) requires a{T) > 0. Otherwise, 
J^{si ^T, fj,) = oo and the simple pole singularity Sg (T, fi) ^ is always the rightmost s-singularity of 2^ ([5]) (see 
Fig. lb). At T > Tc, where (t{T)\smm = 0, the considered system can exist only in the one-phase state. It will be 
shown below that for r > 1 the condition (|14l) can be satisfied even at (j{T) = 0. 

At T < Tc the system undergoes the 1-st order phase transition across the line fi* — ^*{T) defined by Eq. p^ . 
Its explicit form is given by the expression: 



M*(T) 



W{T) 



( TflT 



Z/2 



Tb 



zi exp 



W{T)\ 



T 



k=2 



exp 



a fcV3 
T 



(15) 



-1.5 because it 



The points on the line ^i*{T) correspond to the mixed phase states. First we consider the case r : 
is the standard SMM choice. 

The baryonic density is found as {dp/dfi)T and is given by the following formulae in the liquid and gas phases 



Pi 



dpi 
9/1 



ptd 



respectively. Here the function pid is defined as 



PtdiT,fJ.) = 



/ mT 



3/2 



zi exp 



- bpg 
T 



fci-^ exp 



fe=2 



T 



Due to the condition (|13l) this expression is simplified in the mixed phase: 

3/2 



pTriT) 



mT V 



f 



z\ exp 



k^-^ exp 

k=2 



(16) 



(17) 



(18) 



This formula clearly shows that the bulk (free) energy acts in favor of the composite fragments, but the surface 
term favors single nucleons. 

Since at cr > the sum in Eq. (jl8p converges at any t, is finite and according to Eq. (1161) Pg < 1/b. Therefore, 
the baryonic density has a discontinuity Ap = pi — pg > across the line p*{T) for any r. The discontinuities 
take place also for the energy and entropy densities. The phase diagram of the system in the (T, p)-plane is shown 
in the upper panel of Fig. 2. The line p*{T) (|15p corresponding to the mixed phase states is transformed into the 
finite region in the (T, /9)-plane. As usual, in this mixed phase region of the phase diagram the baryonic density p 
and the energy density are superpositions of the corresponding densities of liquid and gas: 



\ pi + (1 - A) , 



A £; -I- (1 - A) £g 



(19) 



Here A (0 < A < 1) is a fraction of the system volume occupied by the liquid inside the mixed phase, and the 
partial energy densities for (i = l,g) can be found from the thermodynamic identity ^16] : 



.dpH 
dT 



dpi 
dp 



Pi 



(20) 



Inside the mixed phase at constant density p the parameter A has a specific temperature dependence shown in 
the upper panel of Fig. 3: from an approximately constant value p/ Po at small T the function A(r) drops to zero 
in a narrow vicinity of the boundary separating the mixed phase and the pure gaseous phase. This corresponds to 
a fast change of the configurations from the state which is dominated by one infinite liquid fragment to the gaseous 
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multifragment configurations. It liappens inside the mixed phase without discontinuities in the thermodynamical 
functions. 

An abrupt decrease of X{T) near this boundary causes a strong increase of the energy density as a function of 
temperature. This is evident from the middle panel of Fig. 3 which shows the caloric curves at different baryonic 
densities. One can clearly see a leveling of temperature at energies per nucleon between 10 and 20 MeV. As a 
consequence this leads to a sharp peak in the specific heat per nucleon at constant density, Cp(T) = {de/dT)p/p , 
presented in Fig. 3. A finite discontinuity of Cp{T) arises at the boundary between the mixed phase and the gaseous 
phase. This finite discontinuity is caused by the fact that A(T) — 0, but {d\/dT)p 7^ at this boundary (see Fig. 
3). 

It should be emphasized that the energy density is continuous at the boundary of the mixed phase and the 
gaseous phase, hence the sharpness of the peak in Cp is entirely due to the strong temperature dependence of A(T) 
near this boundary. Moreover, at any p < po the maximum value of Cp remains finite and the peak width in Cp(T) 
is nonzero in the thermodynamic limit considered in our study. This is in contradiction with the expectation of 
Refs. |141ll5j that an infinite peak of zero width will appear in Cp(T) in this limit. Also a comment about the 
so-called "boiling point" is appropriate here. This is a discontinuity in the energy density e(T) or, equivalently, a 
plateau in the temperature as a function of the excitation energy. Our analysis shows that this type of behavior 
indeed happens at constant pressure, but not at constant density! This is similar to the usual picture of a liquid-gas 
phase transition. In Refs. [14, 15J a rapid increase of the energy density as a function of temperature at fixed p near 
the boundary of the mixed and gaseous phases (see the lower panel of Fig. 3) was misinterpreted as a manifestation 
of the "boiling point". 

New possibilities appear at non-zero values of the parameter t. At < t < 1 the qualitative picture remains the 
same as discussed above, although there are some quantitative changes. For t > 1 the condition (|T4|) is also satisfied 
at T > Tc where a{T)\sMM = 0. Therefore, the liquid-gas phase transition extends now to all temperatures. Its 
properties are, however, different for r > 2 and for r < 2 (see Fig. 2). If r > 2 the gas density is always lower 
than 1/6 as pid is finite. Therefore, the liquid-gas transition at T > Tc remains the 1-st order phase transition with 
discontinuities of baryonic density, entropy and energy densities (lower panel in Fig. 2) . 

4. The Critical Indices and Scaling Relations of the SMM 

The above results allow one to find the critical exponents a', (3 and 7' of the simplified SMM. These exponents 
describe the temperature dependence of the system near the critical point on the coexistence curve p* — p*{T) 
p^ . where the effective chemical potential vanishes v = p*(T) + W(T) - bp(T,p*(T)) = 



for e < , ^21) 



, for £ > , 

Ap - , for e > , (22) 

KT e"^' , for e > , (23) 

where Ap = pi — pg defines the order parameter, Cp = {'§§') p denotes the specific heat at fixed particle density 
and kt = ^ ^g^^ is the isothermal compressibility. The shape of the critical isotherm for p < Pc is given by the 
critical index 5 (the tilde indicates e — hereafter) 

Pc~P - {pc- for e = . (24) 



The calculation of a and a' requires the specific heat Cp. With the formula [34] 

(25) 



.iT,p) ^ 1 / _ /av 



one obtains the specific heat on the FT curve by replacing the partial derivatives by the total ones [35] . The 
latter can be done for every state inside or on the boundary of the mixed phase region. For the chemical potential 
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/z*(T) = bp*{T) — W{T) one gets ^^^^^ — (p ^ ^) dT^^ ~^ ~~aT^- Here the asterisk indicates the condensation 

hne — 0) hereafter. Fixing p ^ pc = pi = 1/b one finds c*^{T) — T '^ dT^"^ and, hence, obtains a = a' = 0. To 
calculate (3, 7' and S the behavior of the series 

I]q(e,iy) =^ fc*-" eT. " (26) 

fc=2 

should be analyzed for small positive values of e and — // {A = ao/Tc). In the limit e — > the function l]q(£,0) 
remains finite, if t > <? + 1, and diverges otherwise. For t = q + 1 this divergence is logarithmic. The case t < q + 1 
is analyzed in some details, since even in Fisher's papers it was performed incorrectly. 

With the substitution Zk = k [Ae''] one can prove [18] that in the limit £ — ?► the series on the r. h. s. of p6| 
converges to an integral 

00 1 ^ 

l]q(e,0) = {Ae^\~^ ^ -^r^e-^fc {Ae^\~^ f dz z^-^ e~^" . (27) 

k=2 •' 1 

The assumption g — t > — 1 is sufficient to guarantee the convergence of the integral at its lower limit. Using this 
representation, one finds the following general results 

Jir~q~l)^ ifr<g + l, f i^^ ' ^ ' \ if r < g + 1 , 

Sq(e,0)-<| in|£|, ifr = g + l, and Sq(0,;>)~ if r = g + 1 , (28) 

e° , if r>g + l. I if t > g + 1 , 

which allowed us to find out the critical indices of the SMM (see Table 1). 

Table 1. Critical exponents of the SMM and FDM as functions of Fisher index t for the general parametrization 
of the surface energy (T(T)fci e'^fc°' with e = (Tc - T)/Tc 
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In the special case = 2a the well-known exponent inequalities proven for real gases by 



Fisher [36] 
Griffiths[37] 
Liberman[38] 



a' + 2/3 + 7' > 2 , (29) 
a' + /3{l + d) > 2 , (30) 
y + Pil-S) > , (31) 



are fulfilled exactly for any r. (The corresponding exponent inequalities for magnetic systems are often called 
Rushbrooke's, Griffiths' and Widom's inequalities, respectively.) For ( > 2a, Fisher's and Griffiths' exponent 
inequalities are fulfilled as inequalities and for < 2a they are not fulfilled. The contradiction to Fisher's and 
Griffiths' exponent inequalities in this last case is not surprising. This is due to the fact that in the present version 
of the SMM the critical isochore p = Pc = Pi lies on the boundary of the mixed phase to the liquid. Therefore, in 
expression (2.13) in Ref. for the specific heat only the liquid phase contributes and, therefore. Fisher's proof 
of Ref. [36] following (2.13) cannot be applied for the SMM. Thus, the exponent inequalities ((29|) and (pO|) have to 
be modified for the SMM. Using results of Table 141 one finds the following scaling relations 

a' + 2/3 + 7' = - and a' + (3(1 + 6) = - . (32) 
a a 

Liberman's exponent inequality pip is fulfilled exactly for any choice of C and a. 
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Since the coexistence curve of the SMM is not symmetric with respect to p = pc, it is interesting with regard to 
the specific heat to consider the difference Acp{T) = c*^{T) — c*p^{T), following the suggestion of Ref. [35]. Using 
Eq. (P5|) for gas and liquid and noting that l/p* ~b — one obtains a specially defined index a'^ from the most 

divergent term for C > 1 

P^dl^) dT^ [ 2-§((7 + 2-r), if r >(T + 1. 

Then it is a'^ > for ^/cr < 2. Thus, approaching the critical point along any isochore within the mixed phase 
region except for p = pc = 1/b the specific heat diverges for C/f < 2 as defined by a'^ and remains finite for the 
isochore p — Pc = ^/b. This demonstrates the exceptional character of the critical isochore in this model. 

In the special case that C = 1 finds a'^ = 2 — l/cr for r < 1 + 2a and a'^ — —(3 for r > 1 + 2fT. Therefore, 
using a'j, instead of a' , the exponent inequalities ([29)) and ([30]) are fulfilled exactly if C > 1 and t < a + 1 or it ( = 1 
and T < 2(7 + 1. In all other cases and are fulfilled as inequalities. Moreover, it can be shown that the 
SMM belongs to the universality class of real gases for C > 1 and t > a + 1. 

The comparison of the above derived formulae for the critical exponents of the SMM for C = 1 with those 
obtained within the FDM (Eqs. 51-56 in [T9]) shows that these models belong to different universality classes 
(except for the singular case t — 2). 

Furthermore, one has to note that for^ — l,cr<l/2 and 1 + a < t < 1 + 2(7 the critical exponents of the 
SMM coincide with those of the exactly solved one-dimensional FDM with non-zero droplet-volumes |35| . 

For the usual parameterization of the SMM [IJ one obtains with ( — 5/A and a — 2/3 the exponents 

Thus, Fisher suggestion to use a'^ instead of a' allows one to "save" the exponential inequalities, however, it is not 
a final solution of the problem. 

The critical indices of the nuclear liquid-gas FT were determined from the multifragmentation of gold nuclei [39] 
and found to be close to those ones of real gases. The method used to extract the critical exponents /? and 7' in 
Ref. [39] was, however, found to have large uncertainties of about 25 per cents [40]. Nevertheless, those results allow 
us to estimate the value of r from the experimental values of the critical exponents of real gases taken with large 
error bars. Using the above results we generalized [18 the exponent relations of Ref. [3S] 

for arbitrary a and (. Then, one obtains with [IT] P = 0.32 — 0.39 , 7' — 1.3 — 1.4 and 6 = 4 — 5 the estimate 
r = 1.799 — 1.846. This demonstrates also that the value of r is rather insensitive to the special choice of /3 , 7' and 
6, which leads to a'^ = 0.373 — 0.461 for the SMM. Theoretical values for /3 , 7' and S for Ising-like systems within 
the renormalized 4>'^ theory [42] lead to the narrow range t = 1.828 ± 0.001 . The values oi /3 , 7' and S indices for 
nuclear matter and percolation of two- and three-dimensional clusters are reviewed in [23] . 

There was a decent try to study the critical indices of the SMM numerically [44] . The version V2 of Ref. |44) 
corresponds precisely to our model with r = 0, C = 5/4 and <t = 2/3, but their results contradict to our analysis. 
Their results for version V3 of Ref. [44] are in contradiction with our proof presented in Ref. [16]. There it was 
shown that for non- vanishing surface energy (as in version V3) the critical point does not exist at all. The latter 
was found in [33] for the finite system and the critical indices were analyzed. Such a strange result, on one hand, 
indicates that the numerical methods used in Ref. [44] are not self-consistent, and, on the other hand, it shows an 
indispensable value of the analytical calculations, which can be used as a test problem for numerical algorithms. 

It is widely believed that the effective value of r defined as Toff = ~dlnnk{s)/d\nk attains its minimum at 
the critical point (see references in [21]). This has been shown for the version of the FDM with the constraint of 
sufficiently small surface tension a = for T > Tg [43] and also can be seen easily for the SMM. Taking the SMM 
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fragment distribution nk{e) = g{T)k exp[^fc — ^^fc'^] ^ k ''^^ one finds 

TeS = T - ^"'^^ k" =^> r = min(roff) , (36) 

where the last step follows from the fact that the inequalities a{e) > , i' < become equalities at the critical point 
= a(0) — 0. Therefore, the SMM predicts that the minimal value of Tcff corresponds to the critical point, where, 
in contrast to other regions of parameters, the mass distribution of fragments should become power-like. 

In the E900 tt^ + Au multifragmentation experiment [20] the ISiS collaboration measured the dependence of 
Toff upon the excitation energy and found the minimum value min(Teff) = 1.9 (Fig. 5 of Ref. |20| ). Also the 
EOS collaboration [21j performed an analysis of the minimum of Tcff on Au + C multifragmentation data. The 
fitted Toff, plotted in Fig. 11. b of Ref. [21] versus the fragment multiplicity, exhibits a minimum in the range 
min(Tcff) = 1.8 — 1.9 . Both results contradict the original FDM [TO], but agree well with the above estimate of t 
for real gases and for Ising-like systems in general. 

5. Constrained SMM in Finite Volumes 

Despite the great success, the application of the exact solution [T6HT8| to the description of experimental data 
is limited because this solution corresponds to an infinite system and due to that it cannot account for a more 
complicated interaction between nuclear fragments. Therefore, it was necessary to extend the exact solution |16H18) 
to finite volumes. It is clear that for the finite volume extension it is necessary to account for the finite size 
and geometrical shape of the largest fragments, when they are comparable with the system volume. For this one 
has to abandon the arbitrary size of largest fragment and consider the constrained SMM (CSMM) in which the 
largest fragment size is explicitly related to the volume V of the system. Thus, the CSMM assumes a more strict 

K{V) 

constraint ^ k Uk = A , where the size of the largest fragment K{V) = aV/b cannot exceed the total volume of 

k 

the system (the parameter a < 1 is introduced for convenience). The case of the fixed size of the largest fragment, 
i.e. K{V) = Const, analyzed numerically in Ref. |45) is also included in our treatment. A similar restriction should 
be also applied to the upper limit of the product in all partitions Z^(y,T), Za{V,T) and Z{V,T,fi) introduced 
above (how to deal with the real values of K{V), see later). Then the model with this constraint, the CSMM, cannot 
be solved by the Laplace transform method, because the volume integrals cannot be evaluated due to a complicated 
functional F-dependence. However, the CSMM can be solved analytically with the help of the following identity |28] 

GiV) ^ j 1 g e"'(^-«) G{i) , (37) 

— oo — oo 

which is based on the Fourier representation of the Dirac J-function. The representation ([57)) allows us to decouple 
the additional volume dependence and reduce it to the exponential one, which can be dealt by the usual Laplace 
transformation in the following sequence of steps 

+ 00 +00 



i(A,T,^) = dV e-^^ Z{V,T,^i) ^ J^dV J d^ J ^ e^v{v'-i)-xv' 



E 

{"fc} 



fl _L 0fc(r)e '^''^T""^""° }"'' Q{V') = pdV J d^ J ^e^viv'-i)-xv'+v'ni^-^v)_ (gg) 

J^—^ _ —oo — oo 

After changing the integration variable V V ~ V ~ b ^ k Uk, the constraint of 0-function has disappeared. 

k 

Then all rife were summed independently leading to the exponential function. Now the integration over V' in Eq. (|38p 
can be straightforwardly done resulting in 

+ 00 +00 

dr] e-"'« 



Z(A,T,,)^ I I J—^-y^^J—-y (39) 



-OO — oo 
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where the function A) is defined as follows 

K{0 



fe=i 



/ mT 



Zi e 



k^2 



{fi + W-XbT)k- 



(40) 



As usual, in order to find the GCP by the inverse Laplace transformation, it is necessary to study the structure 
of singularities of the isobaric partition 

6. Isobaric Partition Singularities at Finite Volumes 

The isobaric partition PO)) of the CSMM is, of course, more complicated than its SMM analog [16''17] because for 
finite volumes the structure of singularities in the CSMM is much richer than in the SMM, and they match in the 
limit V ^ oo only. To see this let us first make the inverse Laplace transform: 



X+ioc 



z(T/,r,M) 



+ CXD +00 X+*°0 



drj 
2^ 



-00 —00 

+ CXJ +00 

drj 



d\ 



2Tii \ — irj — J-{£,, A — irj) 



d( 



2n 



{A„} 



.A„ V 



dXr, 



(41) 



noindent where the contour A-integral is reduced to the sum over the residues of all singular points A — A^^ 4- ivj 
with n = 1,2,.., since this contour in the complex A-plane obeys the inequality x > max(i?e{A.,i}). Now both 
remaining integrations in (|4ip can be done, and the GCP becomes 

-1 



Z{V,T,t,) 



{A„} 



1 



9A„ 



(42) 



i.e. the double integral in (|1T|) simply reduces to the substitution ^ — > F in the sum over singularities. This is a 
remarkable result which was formulated in Ref. [28' as the following theorem: if the Laplace- Fourier image of the 
excluded volume GCP exists, then for any additional V -dependence of J-{V, Xn) or (pkiT) the GCP can be identically 
represented by Eq. 114 

The simple poles in (|4ip are defined by the equation 



A„ = HV,Xn)- 



(43) 



In contrast to the usual SMM |16l[T7] the singularities A„ are (i) are volume dependent functions, if K{V) is not 
constant, and (ii) they can have a non-zero imaginary part, but in this case there exist pairs of complex conjugate 
roots of because the GCP is real. 

Introducing the real i?„ and imaginary /„ parts of A„ — Rn + iln , we can rewrite Eq. (|43p as a system of coupled 
transcendental equations 

K(V) 

-Rn = E e^^^ cos(/„6A:) , (44) 

k=l 
K{V) 

/„ = - ^ 0fe(r) e^^^^ sin(/„6fc) , (45) 
fc=i 



where we have introduced the set of the effective chemical potentials i/„ = v{\n) with v{\) = /i + W{T) — XhT, 

and the reduced distributions 0i(T) = (^) ^ zi exp(-T^(r)/T) and 4>k>iiT) = (^)^ fc-^ exp(-cr(T) P/^/T) 
for convenience. 

Consider the real root (Rq > 0,Iq =0), first. For /„ = Iq = the real root Rq exists for any T and /i. Comparing 
Rq with the expression for vapor pressure of the analytical SMM solution [T6',T7j shows that TRq is a constrained 
grand canonical pressure of the gas. As usual, for finite volumes the total mechanical pressure [51128] differs from 
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TRq. Equation ([45|) shows that for /n>o 7^ the inequahty cos(/„6A:) < 1 never become the equaUty for all fc-values 
simultaneously. Then from Eq. (|44|) one obtains {n > 0) 

K{V) 



where the second inequality (|46l) immediately follows from the first one. In other words, the gas singularity is 
always the rightmost one. This fact plays a decisive role in the thermodynamic limit V 00. 

The interpretation of the complex roots A„>o is less straightforward. According to Eq. (jH]), the GCP is a 
superposition of the states of different free energies —XnVT. (Strictly speaking, —XnVT has a meaning of the 
change of free energy, but we will use the traditional term for it.) For n > the free energies are complex. Therefore, 
—Xn>oT is the density of free energy. The real part of the free energy density, —RnT, defines the significance of the 
state's contribution to the partition: due to (|46l) the largest contribution always comes from the gaseous state and 
has the smallest real part of free energy density. As usual, the states which do not have the smallest value of the 
(real part of) free energy, i. e. — i?„>oT, are thermodynamically metastable. For infinite volume they should not 
contribute unless they are infinitesimally close to —RqT, but for finite volumes their contribution to the GCP may 
be important. 

As one sees from (|44|) and (|45|). the states of different free energies have different values of the effective chemical 
potential Vn, which is not the case for infinite volume |161I17|. where there exists a single value for the effective 
chemical potential. Thus, for finite V the states which contribute to the GCP (|^^ are not in a true chemical 
equilibrium. 

The meaning of the imaginary part of the free energy density becomes clear from (|44)) and (|45|) |46| : as one can see 
from (|44p the imaginary part In>o effectively changes the number of degrees of freedom of each fc-nucleon fragment 
{k < K{V)) contribution to the free energy density — i?„>oT. It is clear, that the change of the effective number 
of degrees of freedom can occur virtually only and, if Xn>o state is accompanied by some kind of equilibration 
process. Both of these statements become clear, if we recall that the statistical operator in statistical mechanics 
and the quantum mechanical convolution operator are related by the Wick rotation [47J. In other words, the inverse 
temperature can be considered as an imaginary time. Therefore, depending on the sign, the quantity /„6T = 
that appears in the trigonometric functions of the equations (|44|) and (|45p in front of the imaginary time 1/T can 
be regarded as the inverse decay /formation time r„ of the metastable state which corresponds to the pole A„>o (for 
more details see next sections). 

This interpretation of t„ naturally explains the thermodynamic metastability of all states except the gaseous 
one: the metastable states can exist in the system only virtually because of their finite decay/formation time, 
whereas the gaseous state is stable because it has an infinite decay/formation time. 

7. No Phase Transition Case 

It is instructive to treat the effective chemical potential i'(A) as an independent variable instead of /i. In contrast to 
the infinite V, where the upper limit < defines the liquid phase singularity of the isobaric partition and gives 
the pressure of a liquid phase pi{T, /i) = TRo\v^oo = (m + |16[|17| . for finite volumes and finite K{V) the 

effective chemical potential can be complex (with either sign for its real part) and its value defines the number and 
position of the imaginary roots {A„>o} in the complex plane. Positive and negative values of the effective chemical 
potential for finite systems were considered [i 33j within the Fisher droplet model, but, to our knowledge, its complex 
values have never been discussed. From the definition of the effective chemical potential iy{X) it is evident that its 
complex values for finite systems exist only because of the excluded volume interaction, which is not taken into 
account in the Fisher droplet model [19}. However, a recent study of clusters of the d = 2 Ising model within the 
framework of EDM (see the corresponding section in Ref. [23]) shows that the excluded volume correction improves 
essentially the description of the thermodynamic functions. Therefore, the next step is to consider the complex 
values of the effective chemical potential and free energy for the excluded volume correction of the Ising clusters on 
finite lattices. 

As it is seen from the upper panel of Fig. 4, the r.h.s. of Eq. (|45p is the amplitude and frequency modulated 
sine-like function of dimensionless parameter /„ b. Therefore, depending on T and Reiy) values, there may exist 




(46) 
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no complex roots {A„>o}, a finite number of tlrem, or an infinite number of them. In the upper panel of Fig. 4 we 
showed a special case which corresponds to exactly three roots of Eq. (|43| for each value of K{V): the real root 
(Jo = 0) and two complex conjugate roots (±/i). Since the r.h.s. of is monotonously increasing function of 
Re{v), when the former is positive, it is possible to map the T — Re[v) plane into regions of a fixed number of roots 
of Eq. (|43|) . Each curve in the lower panel of Fig. 4 divides the T — Re{v) plane into three parts: for i?e(j^)-values 
below the curve there is only one real root (gaseous phase), for points on the curve there exist three roots, and 
above the curve there are four or more roots of Eq. (|43p . 

For constant values of K{V) = K the number of terms in the r.h.s. of (|45|) does not depend on the volume and, 
consequently, in thermodynamic limit V oo only the rightmost simple pole in the complex A-plane survives out 
of a finite number of simple poles. According to the inequality (|46p . the real root Aq is the rightmost singularity of 
isobaric partition (|39p . However, there is a possibility that the real parts of other roots A,i>o become infinitesimally 
close to i?0 7 when there is an infinite number of terms which contribute to the GCP (|42p . 

Let us show now that even for an infinite number of simple poles in (|42p only the real root Aq survives in the limit 
V — ^ oo. For this purpose consider the limit Re{i'n) 3> T. In this limit the distance between the imaginary parts 
of the nearest roots remains finite even for infinite volume. Indeed, for Reivn) ^ T the leading contribution to the 
r.h.s. of (j45p corresponds to the harmonic with k = K, and, consequently, an exponentially large amplitude of this 
term can be only compensated by a vanishing value of sin (/„ bK), i.e. /„ bK = nn + (S„ with \6n\ <^ tt (hereafter 
we will analyze only the branch /„ > 0), and, therefore, the corresponding decay/formation time t„ ~ K[TmT]^^ 
is volume independent. 

Keeping the leading term on the r.h.s. of (|45|) and solving for 5m one finds 



/„ « (-ir^0K(T)e^J„, with ^.«^-S^e-^, (47) 

Kb (pK(T) 

Rn « (-l)"^i^(T) e^^^^ , (48) 

where in the last step we used Eq. (j44|) and condition |5„| ^ tt. Since for F ^> oo all negative values of i?„ cannot 
contribute to the GCP (|^^ . it is sufficient to analyze even values of n which, according to (HH), generate i?„ > 0. 

Since the inequality ()46p can not be broken, a single possibility, when A„>o pole can contribute to the partition 
(|42|) . corresponds to the case _R„ Rq — for some finite n. Assuming this, we find i?e(z^(A„)) — > i?e(i/(Ao)) for 
the same value of /i. Substituting these results into equation (H^ . one gets 



K 

Rn w 2^0fc(T) e T cos 

k=l 



nnk 



K 



< Ro . (49) 



The inequality follows from the equation for i?o and the fact that, even for equal leading terms in the sums 
above (with k = K and even n), the difference between Rq and i?„ is large due to the next to leading term k = K — 1, 

-Re(i^(AQ)) (Jf-1) . ... . I 

which is proportional to e t ^1. Thus, we arrive at a contradiction with our assumption Rq — Rn — > 0+, 

and, consequently, it cannot be true. Therefore, for large volumes the real root Aq always gives the main contribution 
to the GCP (|^^ . and this is the only root that survives in the limit V — oo. Thus, we showed that the model with 
the fixed size of the largest fragment has no phase transition because there is a single singularity of the isobaric 
partition (|39p . which exists in thermodynamic limit. 

8. Finite Volume Analogs of Phases 

If K{V) monotonically grows with the volume, the situation is different. In this case for positive value of Re{v) ^ T 
the leading exponent in the r.h.s. of (|45p also corresponds to a largest fragment, i.e. to fc = KiV). Therefore, we can 
apply the same arguments which were used above for the case K{V) = K = canst and derive similarly equations 
((Tf|) - (P5)) for /„ and _R„. From w ^ k(v) follows that, when V increases, the number of simple poles in (|^T|) 
also increases and the imaginary part of the closest to the real A-axis poles becomes very small, i.e /„ — > for 
n ^ K{V), and, consequently, the associated decay/formation time t„ « K{V)[KnT]~'^ grows with the volume of 
the system. Due to /„ — ?• 0, the inequality (|49| cannot be established for the poles with n <^ K{V). Therefore, in 
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contrast to the previous case, for large K{V) the simple poles with n ^ K{V) will be infinitesimally close to the 
real axis of the complex A-plane. 
From Eq. (gH]) it follows that 



T 



1 



K{V)h 



In 



Rn 



Pi [T, fJ-) 
T 



(50) 



for \fi\ ^ T and K{V) oo. Thus, we proved that for infinite volume the infinite number of simple poles moves 
toward the real A-axis to the vicinity of liquid phase singularity A; = pi{T, fJ.)/T of the isobaric partition |161I17| and 
generates an essential singularity of function J-{V,pi/T) in (j40p irrespective to the sign of the chemical potential /i. 
In addition, as we showed above, the states with Re{v) ^ T become stable because they acquire infinitely large 
decay /formation time r„ in the limit V ^ oo. Therefore, these states should be identified as a liquid phase for finite 
volumes as well. 

Now it is clear that each curve in the lower panel of Fig. 4 is the finite volume analog of the phase boundary 
T — /i for a given value of K{V): below the phase boundary there exists a gaseous phase, but at and above each curve 
there are states which can be identified with a finite volume analog of the mixed phase, and, finally, at Re{v) S> T 
there exists a liquid phase. When there is no phase transition, i.e. K{V) ~ K ^ const, the structure of simple poles 
is similar, but, first, the line which separates the gaseous states from the metastable states does not change with 
the volume, and, second, as shown above, the metastable states will never become stable. Therefore, a systematic 
study of the volume dependence of free energy (or pressure for very large V) along with the formation and decay 
times may be of a crucial importance for experimental studies of the nuclear liquid gas phase transition. 

The above results demonstrate that, in contrast to Hill's expectations |[6j, the finite volume analog of the mixed 
phase does not consist just of two pure phases. The mixed phase for finite volumes consists of a stable gaseous phase 
and the set of metastable states which differ by the free energy. Moreover, the difference between the free energies 
of these states is not surface- like, as Hill assumed in his treatment [6], but volume- like. Furthermore, according to 
Eqs. ([33]) and (|45p . each of these states consists of the same fragments, but with different weights. As seen above for 
the case Re{i>) ^ T, some fragments that belong to the states, in which the largest fragment is dominant, may, in 
principle, have negative weights (effective number of degrees of freedom) in the expression for i?„>o (|44p . This can 
be understood easily because higher concentrations of large fragments can be achieved at the expense of the smaller 
fragments and is reflected in the corresponding change of the real part of the free energy —Rn>oVT. Therefore, the 
actual structure of the mixed phase at finite volumes is more complicated than was expected in earlier works. 

The Hills' ideas were developed further in Ref. [5], where the authors claimed to establish the one to one 
correspondence between the bimodal structure of the partition of measurable quantity B known on average and 
the properties of the Lee- Yang zeros of this partition in the complex g-plane. The starting point of Ref. [5j is to 
postulate the partition Zg and the probability Pg{B) of the following form 



Zn = I dB W{B) e-^-s 



Pg{B) ^ 



W{B) e-^ a 



Z„ 



(51) 



where W(B) is the partition sum of the ensemble of fixed values of the observable {i?} , and g is the corresponding 
Lagrange multiplier. Then the authors of Ref. [5j assume the existence of two maxima of the probability Pg{B) (= 
bimodality) and discuss their relation to the Lee- Yang zeros of Zg in the complex g-plane. 

The CSMM gives us a unique opportunity to verify the Chomaz and Gulminelli idea on the bimodality behavior 
of Pg{B) using the first principle results. Let us use the equation ([55]) identifying the intensive variable g with A 



and extensive one B with the available volume V' V. The evaluation of the r.h.s. of (|38l) is very difficult in 
general, but for a special case, when the eigen volume b is small this can be done analytically. Thus, approximating 
J'(^, X-iri) K A) - i7] dT{^, X)/dX, we obtain the CSMM analog of the probability ^ 



-\-oo 



P^{V) Z{\,T,ii)= / 



27r 



(52) 
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where we made the rj integration after applying the approximation for J-{£,., X—irj). Further evaluation of ([52|) requires 
to know all possible solutions of the average volume of the system £,a{V) = V — , X)/dX {a = {1, 2, . . . }). It 

can be shown [46J that for the gaseous domain ly = Re{v) < —2T (see the lower panel of Fig. 4) there exist a single 
solution a — 1, whereas for the domain — Re{v) > 0, which corresponds to a finite volume analog of the mixed 
phase, there are two solutions a = 1, 2. In contrast to the expectations of Ref. [5], the probability ((52)) 



1 ,y[^(Ca)-A] ^ dlnP.jV) 



1 



dV 



< 0, (53) 



has negative derivative for the whole domain of existence of the isobaric partition Z{X, T, ii) [46J. This is true even for 
the domain in which, as we proved, there exists a finite analog of the mixed phase, i.e. for ;/ — Reiy) > 0. Moreover, 

irrespective to the sign of the derivative ■, the probability cannot be measured experimentally. Above 

it was rigorously proven that for any real f the IP Z(A, T, /i) is defined on the real A-axis only for J^(f , A) — A > 0, 
i.e. on the right hand side of the gaseous singularity Aq: A > Aq. However, as one can see from the equation (HI]), 
the "experimental" A„ values belong to the other region of the complex A-plane: i?e(A„>o) < Aq. 

Thus, it turns out that the suggestion of Ref. to analyze the probability (|5ip does not make any sense 
because, as we showed explicitly for the CSMM, it cannot be measured. It seems that the starting point of the 
Ref. [S] approach, i.e. the assumption that the left equation (j5ip gives the most general form of the partition 
of finite system, is problematic. Indeed, comparing ([50)) with the analytical result (|55)) . we see that for finite 
systems, in contrast to the major assumption of Ref. [5j, the probability W of the CSMM depends not only on the 
extensive variable V ^ but also on the intensive variable A, which makes unmeasurable the whole construct of Ref. [5]. 
Consequently, the conclusions of Ref. [5l on the relation between the bimodality and the phase transition existence 
are not general and have a limited range of validity. In addition, the absence of two maxima of the probability ()53p 
automatically means the absence of back-banding of the equation of state [5J . 

9. Gas of Bags in Finite Volumes 

Now we will apply the formalism of the preceding sections to the analysis of the Gas of Bags Model (GBM) |31lH5] 
in finite volumes. In the high and low temperature domains the GBM reduces to two well known and successful 
models: the hadron gas model |j49j and the bag model of QGP [50J. Both of these models are surprisingly successful 
in describing the bulk properties of hadron production in high energy nuclear collisions, and, therefore, one may 
hope that their generalization, the GBM, may reflect basic features of the nature in the phase transition region. 
The van der Waals gas consisting of n hadronic species, which are called bags in what follows, has the following 

GCP 131] 

^iy-T)=Y.\i\ i<'^--'^'--y-'^")^'<^)i"1 ..-.„«,, , (54) 

where 4>k{T) = 4>{T,mk) = ^ P^dp exp [- {p^ + mlf'^ /T] = K2 (^) is the particle density 

of bags of mass mfc and eigen volume Vk and degeneracy gk ■ This expression differs slightly form the GCP of the 
simplified SMM ([3]) , where = and the eigen volume of /c-nucleon fragment kb is changed to the eigen volume of 
the bag Vk ■ Therefore, as for the simplified SMM the Laplace transformation ^ with respect to volume of Eq. (|5^ 
gives 



Z(s,T) 



(55) 



In preceding sections we showed that as long as the number of bags, n, is finite, the only possible singularities of 
Z[s, T) ()55[) are simple poles. However, in the case of an infinite number of bags an essential singularity of Z(s, T) 
may appear. This property is used the GBM: the sum over different bag states in ([54)) can be replaced by the 
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integral, X^j^iffi--- — dmdv...p{m,v), if the bag mass-volume spectrum, p{rn,v), which defines the number 
of bag states in the mass- volume region [to, t;;m -f dm,v + dv], is given. Then, the Laplace transform of Z{V,T) 
reads ^31] 

oo oo 

Zgb{s,T) = J dV e""^ Z(y,T) = [s - /(r,s)]"^ , where /(r,s) = J dmdv p(m,w) e"'"' </>(T, m) . (56) 



Like in the simplified SMM, the pressure of infinite system is again given by the rightmost singularity: p{T) = 
Ts*{T) = T ■ max{sH{T), sq{T)}. Similarly to the simplified SMM considered in Sect. II and III, the rightmost 
singularity s*{T) of Z{s,T) can be either the simple pole singularity sh{T) — f {T, sh{T)) of the isobaric 
partition ([TO]) or the sq{T) singularity of the function f{T,s) ([55]) it-self. 

The major mathematical difference between the simplified SMM and the GBM is that the latter employs the 
two parameters mass-volume spectrum. Thus, the mass-volume spectrum of the GBM consists of the discrete 
mass- volume spectrum of light hadrons and the continuum contribution of heavy resonances j51| 

p(to,w)=^^ gj S{rn~~'mj) S{v — Vj) + Q{v — Vo)Q{'m~ AIo~~ Bv)C v'^ {m — Bv)^ exp 
i=i 

respectively. Here nij < Mq, Vj < Vq, Mq « 2 GeV, Vp ~ 1 fm'^, C, 7, S and B (the so-called bag constant, B « 400 
MeV/fm'^) are the model parameters and 

TT^ / 7 \ TT^ / 7 \ TT^ 95 

'^Q = ^ 9, + ^9,,] = ^h-S + -•2.2.3.3 = - - (58) 



4 11 ; 

— ag (m — Bv)i 
3 



(57) 



30 V 8"'7 30 V 8 y 30 2 

is the Stefan-Boltzmann constant counting gluons (spin, color) and (anti-)quarks (spin, color and u, d, s-flavor) 
degrees of freedom. 

Recently the grand canonical ensemble has been heavily criticized [SH[53], when it is used for the exponential 
mass spectrum. This critique, however, cannot be applied to the mass-volume spectrum (|57l) because it grows 
less fast than the Hagedorn mass spectrum discussed in [52l[53j and because in the GBM there is an additional 
suppression of large and heavy bags due to the van der Waals repulsion. Therefore, the spectrum ([ST)) can be safely 
used in the grand canonical ensemble. 

It can be shown |48| that the spectrum ((57|) generates the sq{T) = ~ singularity, which reproduces 

the bag model pressure p{T) = Tsq{T) [50] for high temperature phase, and sh{T) singularity, which gives the 
pressure of the hadron gas model [IS] for low temperature phase. The transition between them can be of the first 
order or second order or cross-over, depending on the model parameters. 

However, for finite systems the volume of bags and their masses should be finite. The simplest finite volume 
modification of the GBM is to introduce the volume dependent size of the largest bag n = n{V) in the partition 
(|54p . As we discussed earlier such a modification cannot be handled by the traditional Laplace transform technique 
used in , but this modification can be easily accounted for by the Laplace-Fourier method ^28j . Repeating all 

the steps of the sections V and VI, we shall obtain the equations (|40|) - (|43)) . in which the function J-'(^, A) should be 
replaced by its GBM analog /(A, Vb) = ./ff(A) + /q(A, Vb) defined via 

,/,„ Vb/Vo 

fniX) = 9j <t>{T,m,) e""^^ , and /q(A, Vb) = Vo j dk a{T, Vok) e^°(^«(^)-^)'= . (59) 
3=1 i 

In evaluating (j59p we used the mass- volume spectrum (|57p with the maximal volume of the bag Vb and changed 
integration to a dimensionless variable k — v/Vq. Here the function a{T,v) = m(T)u^+''+'' is defined by u{T) — 
C^-i^+'Z' T^+^^iagT^ + Bf/^. 

The above representation (j59p generates equations for the real and imaginary parts of A„ = i?„ + i/„, which 
are very similar to the corresponding expressions of the CSMM ([^^ and (|^5l) . Comparing ([5^)1 with one sees 
that their main difference is that the sum over k in ((43|) is replaced by the integral over k in (|59p . Therefore, the 
equations (j44| and (|45l) remain valid for i?„ and In of the GBM, respectively, if we replace the k sum by the 
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integral for K{V) = Vb/Vo, b = Vo, t^(A) = Vo{sq{T) - A) and <^fe>i(r) = Vo a{T,Vok). Thus, the results and 
conclusions of our analysis of the i?„ and /„ properties of the CSMM should be valid for the GBM as well. In 
particular, for large values of K{V) = Vb/Vq and i?„ < sq{T) one can immediately find out /„ w -kti/Vb and the 
GBM formation/decay time t„ = Vb [TrnTVb]^^. These equations show that the metastable A„>o states can become 
stable in thermodynamic limit, if and only if Vb ^ V . 

The finite volume modification of the GBM equation of state should be used for the quantities which have 
V\o ^ 1. This may be important for the early stage of the relativistic nuclear coUisons when the volume of the 
system is small, or for the systems that have small pressures. The latter can be the case for the pressure of strange 
or charm hadrons. 

10. Hills and Dales Model and the Source of Surface Entropy 

During last forty years the Fisher droplet model (FDM) [12] has been successfully used to analyze the condensation 
of a gaseous phase (droplets or clusters of all sizes) into a liquid. The systems analyzed with the FDM are many 
and varied, but up to now the source of the surface entropy is not absolutely clear. In his original work Fisher 
postulated that the surface free-energy F4 of a cluster of A-constituents consists of surface [A^^^) and logarithmic 
(InA) parts, i.e. Fa = a{T) A^/^ + tTIuA. Its surface part a{T) ^^/s = ^^[^ _ y/y^j ^2/3 consists of the surface 
energy, i.e. ctq A^^^, and surface entropy —Oo/Tc A^/^ . From the study of the combinatorics of lattice gas clusters 
in two dimensions. Fisher postulated the specific temperature dependence of the surface tension <7(T)|fdm which 
gives naturally an estimate for the critical temperature T^. Surprisingly Fisher's estimate works for the 3-d Ising 
model [54], nucleation of real fiuids [55l[56], percolation clusters [57] and nuclear multifragmentation [3]. 

To understand why the surface entropy has such a form we formulated a statistical model of surface deformations 
of the cluster of A-constituents, the Hills and Dales Model (HDM) [29]. For simplicity we consider cylindrical 
deformations of positive height hk > (hills) and negative height —hk (dales), with fc-constituents at the base. It 
is assumed that cylindrical deformations of positive height hk > (hills) and negative height —h^ (dales), with 
fc-constituents at the base, and the top (bottom) of the hill (dale) has the same shape as the surface of the original 
cluster of ^-constituents. We also assume that: (i) the statistical weight of deformations exp (— cro|A5'fc|/si/r) is 
given by the Boltzmann factor due to the change of the surface |A5fc| in units of the surface per constituent si; 
(ii) all hills of heights hk < Hk {Hk is the maximal height of a hill with a base of fc-constituents) have the same 
probability dhk/Hk besides the statistical one; (iii) assumptions (i) and (ii) are valid for the dales. 

The HDM grand canonical surface partition (GCSP) 



Z{Sa 



Kr, 



n 

fc=i 



e(sig) 



(60) 



corresponds to the conserved (on average) volume of the cluster because the probabilities of hill z'^ and dale z^. of 



the same fc-constituent base are identical [29 



dhk -L 



±H^ 



k 



(JoPkHk 



1 - e 



(61) 







Here Pk is the perimeter of the cylinder base. 

The geometrical partition (degeneracy factor) of the HDM or the number of ways to place the center of a 



given deformation on the surface of the j4-constituent cluster which is occupied by the set of {r 
deformations of the ^-constituent base we assume to be given in the van der Waals approximation 



Q = 



K„ 



E k{nl 



k=l 



1^1 



0,1,2,...} 



(62) 



where sik is the area occupied by the deformation of /c-constituent base (fc = 1, 2, ...), Sa is the full surface of the 
cluster, and K^axiSA) is the A-dependent size of the maximal allowed base on the cluster. 

The 0(siC/)-function in (|T|) ensures that only configurations with positive value of the free surface of cluster 
are taken into account, but makes the analytical evaluation of the GCSP (|T]) very difficult. However, we were able 
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to solve this GCSP exactly for any surface dependence of Kmax{SA) using identity ([37]) of the Laplace- Fourier 
transform technique [28]: 

Z{Sa) = - ^^^fe^l ■ (63) 



The poles A„ of the isochoric partition are defined by 

g-fc.iA„^ (64) 



A,i = J"(5a,A„)= ^ 
fc=i 



Sl Si 



Our analysis shows that Eq. (O has exactly one real root i?o = Aq, /to(Ao) = 0, which is the rightmost singularity, 
i.e. Rq > i?e(A„>o). As proved in [29J, the real root Rq dominates completely for clusters with A > 10. 

Also we showed that there is an absolute supremum for the real root Rq , which corresponds to the limit of infi- 
nitesimally small amplitudes of deformations, Hk — >■ 0, of large clusters: sup(i?o) — 1-06009 = Rq = 2 [e^" — l] 
It is remarkable that the last result is, first, model independent because in the limit of vanishing amplitude of 
deformations all model specific parameters vanish; and, second, it is valid for any self-non-intersecting surfaces. 

For large spherical clusters the GCSP becomes Z{Sa) ~ 0.3814 el-06009A ^ which, combined with the 
Boltzmann factor of the surface energy e~'^° generates the following temperature dependent surface tensi- 



on of the large cluster <j{T) — ctq 



1 - 1.06009^ 



This result means that the actual critical temperature of the 



FDM should be Tc = (Tq/I -06009, i.e. 6.009 % smaller in CTq units than Fisher originally supposed. 
11. Strategy of Success 

Here we discussed exact analytical solutions of a variety of statistical model which are obtained by a new powerful 
mathematical method, the Laplace-Fourier transform. Using this method we solved the constrained SMM and Gas 
of Bags Model for finite volumes, and found the surface partition of large clusters. Since in the thermodynamic limit 
the CSMM has the nuclear liquid-gas PT and the GEM describes the PT between the hadron gas and QGP, it was 
interesting and important to study them for finite volumes. As we showed, for finite volumes their GCP functions 
can be identically rewritten in terms of the simple poles A„>o of the isobaric partition (|39p . We proved that the real 
pole Ao exists always and the quantity TAq is the constrained grand canonical pressure of the gaseous phase. The 
complex roots A„>o appear as pairs of complex conjugate solutions of equation (|43p . Their most straightforward 
interpretation is as follows: — Ti?e(A„) has a meaning of the free energy density, whereas bTIm{Xn), depending on 
its sign, gives the inverse decay/formation time of such a state. Therefore, the gaseous state is always stable because 
its decay/formation time is infinite and because it has the smallest value of free energy, whereas the complex poles 
describe the metastable states for i?e(A„>o) > and mechanically unstable states for i?e(A„>o) < 0. 

We studied the volume dependence of the simple poles and found a dramatic difference in their behavior for the 
case with phase transition and without it. For the case with phase transition this formalism allows one to define the 
finite volume analogs of phases unambiguously and to establish the finite volume analog of the T — ^ phase diagram 
(see Fig. 4). At finite volumes the gaseous phase is described by a simple pole Aq, the mixed phase corresponds to 
a finite number of simple poles (three and more) , whereas the liquid is represented by an infinite amount of simple 
poles at -> oo which describe the states of a highest possible particle density. 

As we showed for the CSMM and GEM, at finite volumes the A„ states of the same partition with given T 
and /i are not in a true chemical equilibrium because the interaction between the constituents generates complex 
values of the effective chemical potential. This feature cannot be obtained within the Fisher droplet model due to 
lack of the hard core repulsion between the constituents. We showed that, in contrast to Hill's expectations [6J, the 
mixed phase at finite volumes is not just a composition of two states which are the pure phases. As we showed, a 
finite volume analog of the mixed phase is a superposition of three and more collective states, and each of them is 
characterized by its own value of A„ , and, consequently, the difference between the free energies of these states is 
not a surface-like, as Hill argued but volume-like. 
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Also the exact analytical formulas gave us a unique opportunity to verify the Chomaz and Gulminelli ideas [5] 
about the connection between bimodality and the phase transition existence for finite volumes. The CSMM exact 
analytical solution not only provided us with a counterexample for which there is no bimodality in case of finite 
volume phase transition, but it gave us an explicit example to illustrate that the probability which, according to 
Ref. [5] is supposed to signal the bimodal behavior of the system, cannot be measured experimentally. 

All this clearly demonstrates that the exactly solvable models are very useful theoretical tools and they open 
the new possibilities to study the critical phenomena at finite volumes rigorously. The short range perspectives 
(SRP) of this direction of research are evident: 

1. Study the isobaric ensemble and the excluded volume correction for the clusters of the 2- and 3-dimensional 
Ising models, and find out the reliable signals of phase transition on finite lattices. 

2. Widen or refine the CSMM and GMB analytical solutions for more realistic interaction between the consti- 
tuents. In particular, a more realistic Coulomb interaction between nuclear fragments (not the Wigner- 
Seitz one!) can be readily included now into the CSMM and may be studied rigorously without taking 
thermodynamic limit. 

3. Deepen or extend the CSMM and GMB models to the canonical and microcanonical formulations, and 
work out the reliable signals of the finite system phase transitions for this class of models. 

The major goals for the SRP are (I) to get the reliable experimental signals obtained not with the ad hoc 
theoretical constructs which are very popular nowadays, but directly from the first principles of statistical mechanics; 
(II) to work out a common and useful theoretical language for a few nuclear physics communities. 

However, even the present (very limited!) amount of exact results can be used as a good starting point to build 
up a truly microscopic theory of phase transitions in finite systems. In fact, the exact analytical solution, which we 
found for finite volumes, is one of the key elements that are necessary to create a microscopic kinetics of PTs in 
finite systems. The formulation of such a theory for nuclear physics is demanded by the reality of the experimental 
measurements: both of the phase transitions which are studied in nuclear laboratories, the liquid-gas and hadron 
gas - QGP, are accessible only via the violent nuclear collisions. As a result, in these collisions we are dealing with 
the PTs which occur not only in finite system, but in addition these PTs happen dynamically. It is known that 
during the course of collision the system experiences a complicated evolution from a highly excited (on the ordinary 
level) state which is far from local equilibrium, to the collective expansion of the locally thermalized state and to a 
(nearly) free-streaming stage of corresponding constituents. 

A tremendous complexity of the nuclear collision process makes it extremely difficult for theoretical modeling. 
This is, in part, one of the reasons why, despite a great amount of experimental data collected during last 25 years 
and numerous theoretical attempts, neither the liquid-gas nor the hadron gas - QGP phase transitions are well 
established experimentally and well understood theoretically. It turns out that the major problem of modeling both 
of these PTs in dynamic situations is the absence of the suitable theoretical apparatus. 

For example, it is widely believed in the Relativistic Heavy Ion community (RHIc) that relativistic hydrodynami- 
cs is the best theoretical tool to model the PT between QGP and hadron gas because it employs only the equilibrium 
equation of state [SH]. Up to now this is just a wishful illusion because besides the incorrect boundary conditions, 
known as freeze-out procedure" |591 l60j. which are typically used in the actual hydro calculations [55], the employed 
equation of state does not fit into the finite (and sometimes small!) size of the system because it corresponds to an 
infinite system. On the other hand it is known |61| that hydrodynamic description is limited by the weak (small) 
gradients of the hydro variables, which define a characteristic scale not only for collective hydro properties, but also 
a typical volume for the equation of state. 

Above we showed that for finite systems the equation of state inevitably includes the volume dependence of such 
thermodynamic variables as pressure and energy density which are directly involved into hydrodynamic equations. 
This simple fact is not realized yet in the RHIc, but, probably, the chemical non-equilibrium (which is usually 
implemented into equation of state by hand) is, in part, generated by the finite volume corrections of the GCP. If 
this is the case, then, according to our analysis of the finite volume GCP functions, it is necessary to insert the 
complex values of the chemical potential into hydro calculations. 
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Unfortunately, at present there is no safe recipe on how to include the finite volume equation of state in the 
hydrodynamic description. A partial success of the hybrid hydro-cascade models |62[I63|. which might be considered 
as a good alternative to hydrodynamics, is compensated by the fact that none of the existing hydro-cascade models 
was able to resolve the so called HBT puzzles [13] found in the energy range of the Relativistic Heavy Ion Collider. 
Moreover, despite the rigorous derivation [641165] of the hydro-cascade equations, the hydrodynamic part of this 
approach is suffering the very same problems of the infinite matter equation of state which we discussed above. 
Therefore, further refinements of the hydro-cascade models will not be able to lift up the theoretical apparatus 
of modeling the dynamics of the finite volume PTs to new heights, and we have to search for a more elaborate 
approach. 

It turns out that the recently derived finite domain kinetic equations |64ll65j can provide us with another starting 
point to develop the first principle microscopic theory of the critical phenomena in finite systems. These equations 
generalize the relativistic Boltzmann equation to finite domains and, on one hand, allow one to conjugate two 
(different!) kinetics which exist in two domains separated by the evolving boundary, and, on the other hand, to 
account exactly for the exchange of particles between these domains. (For instance, one can easily imagine the 
situation when on one side of the boundary separating the domains there may exist one phase of the system which 
interacts with the other phase located on the other side of the boundary.) But, first, the finite domain kinetic 
equations should be generalized to the two-particle distribution functions and then they should be adapted to the 
framework of nuclear multifragmentation and the Gas of Bags Model. In doing this, the exact analytical results we 
discussed will be indispensable because they provide us with the equilibrium state of the finite system and tell us 
to what finite volume analog of phases this state belongs. 

Therefore, a future success in building up a microscopic kinetics of PTs in finite systems can be achieved, if 
we combine the exact results obtained for equilibrated finite systems with the rigorous kinetic equations suited for 
finite systems. There is a good chance for the nuclear multifragmentation community to play a very special role in 
the development of such a theory, namely it may act as a perfect and reliable test site to work out and verify the 
whole concept. This is so because besides some theoretical advances and experience in studying the PTs in finite 
systems, the experiments at intermediate energies, compared to the searches for QGP, are easier and cheaper to 
perform, and the PT signals are cleaner and unspoiled by a strong flow. Moreover, once the concept is developed 
and verified, it can be modified and applied to study other PTs in finite systems, including the transitions to/from 
high temperature QCD and dense hadronic matter planned to be studied at CERN LHC and GSI FAIR. Thus, 
after some readjustment the manpower and experimental facilities of nuclear multifragmentation community can 
be used for a new strategic aim, which is at the frontier line of modern physics. 

The fruitful discussions with J. B. Elliott, L.G. Moretto and L. Phair are appreciated. This work was in part 
supported by the US Department of Energy. P.T.R. acknowledges the financial support from the German Academic 
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Fig. 2. Phase diagram in T — p plane for r = —1.5 (upper panel), r = 1.1 (middle panel) and r = 2.1 (lower panel). The mixed phase is 
represented by the extended region. Liquid phase (shown by crosses) exists at density p = po- Point C in the upper panel is the critical 
point, whereas in the middle panel it is the tricritical point. For r > 2 (lower panel) the PT exists for all temperatures T > 0. 
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Fig. 3. Upper panel: Volume fraction A(T) of the liquid inside the mixed phase is shown as a function of temperature for fixed nucleon 
densities p/po = 1/6, 1/3, 1/2, 2/3, 5/6 (from bottom to top) and r = —1.5. 

Middle panel: Temperature as a function of energy density per nucleon (caloric curve) is shown for fixed nucleon densities p/po = 
1/6, 1/3, 1/2, 2/3 and r = —1.5. Note that the shape of the model caloric curves is very similar to the experimental finding 32 , although 
our estimates for the excitation energy is somewhat larger due to oversimplified interaction. For a quantitative comparison between the 
simplified SMM the full SMM interaction should be accounted for. 

Lower panel: Specific heat per nucleon as a function of temperature at fixed nucleon density p/po = 1/3. The dashed line shows the 
finite discontinuity of Cp(T) at the boundary of the mixed and gaseous phases for r = —1.5. 
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Fig. 4. Upper panel: A graphical solution of Eq. H45| l for T = 10 MeV and r = 1.825. The l.h.s. (straight line) and r.h.s. of Eq. ||45|I 
(all dashed curves) are shown as the function of dimensionless parameter Ii b for the three values of the largest fragment size K(V). 
The intersection point at (0; 0) corresponds to a real root of Eq. 1431 1. Each tangent point with the straight line generates two complex 
roots of l [43t . 

Lower panel: Each curve separates the T — Reivn) region of one real root of Eq. 1)43^ (below the curve), three complex roots (at the 
curve) and four and more roots (above the curve) for three values of K{V) and the same parameters as in the upper panel. 
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